Implementation of low-loss superinductances for quantum circuits 
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The simultaneous suppression of charge fluctuations and offsets is crucial for preserving quantum 
coherence in devices exploiting large quantum fluctuations of the superconducting phase. This re- 
quires an environment with both extremely low DC and high RF impedance. Such an environment is 
provided by a superinductance, defined as a zero DC resistance inductance whose impedance exceeds 
the resistance quantum Rq — h/(2e) 2 ~ 6.5 kfl at frequencies of interest (1 - 10 GHz). In addition, 
the superinductance must have as little dissipation as possible, and possess a self-resonant frequency 
well above frequencies of interest. The kinetic inductance of an array of Josephson junctions is an 
ideal candidate to implement the superinductance provided its phase slip rate is sufficiently low. We 
successfully implemented such an array using large Josephson junctions (Ej » Ec), and measured 
internal losses less than 20 ppm, self-resonant frequencies greater than 10 GHz, and phase slip rates 
less than 1 mHz. 

PACS numbers: 85.25.Cp, 74.81. Fa, 74.50.+r, 64.70.Tg 



The emerging field of quantum electronics exploiting 
large fluctuations of the superconducting phase is lim- 
ited by the practical challenge of engineering an electro- 
magnetic environment which suppresses simultaneously 
the quantum fluctuations of charge and the random low- 
frequency fluctuations of offset charges. The small value 
of the fine structure constant a — 1/137 entails a fun- 
damental asymmetry between flux and charge quantum 
fluctuations, strongly favoring the latter. To illustrate 
this, let us consider the simplest case of a dissipationless 
LC oscillator, where the charge Q on the capacitor plates 
and the generalized flux $ across the inductor are conju- 
gate variables. The ratio between quantum fluctuations 
of charge, 5q = 8Q/(2e), and flux, Sip = (5$/$ , in the 
ground state of the oscillator is given by Sip/Sq = Z /Rq. 
Here Z = \jLjC is the characteristic impedance of the 
oscillator and Rq = h/(2e) 2 = 6.5fcf2 is the supercon- 
ducting resistance quantum. Using only geometrical in- 
ductors and capacitors, the characteristic impedance of 
the oscillator Zq can not exceed the vacuum impedance 
Z V ac — V^oAoj thus imposing quantum fluctuations of 
charge at least an order of magnitude larger than flux 
fluctuations: 8ip/8q < Z vac /RQ = 8a. 

In order to exceed the vacuum impedance and sup- 
press quantum charge fluctuations, circuit elements with 
extremely high impedance are required. On chip resistors 
and long chains of Josephson junctions (JJs) in the dis- 
sipative regime have already been used to provide high 
impedance environments [THH] - However, these Ohmic 
components do not help shunting charge offsets and, be- 
ing dissipative, they tend to destroy the quantum coher- 
ence of the devices. In order to simultaneously suppress 
quantum charge fluctuations and offsets, we need a cir- 
cuit clement which possesses three key attributes: high 
impedance at frequencies of interest, perfect conduction 
at DC and extremely low dissipation. These attributes 
define the so-called "superinductance" [I]. 

In this Letter, we report the successful implementa- 



tion and detailed characterization of superinductances 
using the large kinetic inductance of arrays of Josephson 
junctions (see Fig[lja)), following the proof-of-concept 
shown in the fluxonium circuit [5] . This first superinduc- 
tance implementation suffered from coherent quantum 
phases-slips (CQPS) [5] [7J, which constitute additional 
degrees of freedom, difficult to control experimentally. 
We show that we can completely suppress the CQPS 
by using large JJs with Josephson energy Ej^ 100 Ec, 
where Ec — e 2 / (2Cj) is the charging energy of one junc- 
tion. 

Superconducting nanowires are also likely candidates 
for implementing superinductances [5J [5] • Unfortunately, 
they show significant internal dissipation which is not yet 
well understood, and are more challenging to fabricate. 
The JJ arrays may also suffer from dissipation, either 
due to coupling to internal degrees of freedom [lOl [11] 
or to a dissipative external bath [12]. Indeed, transport 
measurements on large arrays of JJ show the appearance 
of a superconducting to insulating transition (SIT) with 
decreasing Josephson energy Ej [PSl - fTS] . Previous im- 
plementations of resonators where the inductive energy 
is given by JJ arrays have yielded internal quality factors 
in the range of a few thousands [151 EZ] • We show here 
that the problem of dissipation can be alleviated by using 
large JJs to realize superinductances far from the SIT re- 
gion with internal quality factors an order of magnitude 
above the previously reported values. 

The superinductances are formed by an array of closely 
spaced Josephson junctions on a C-plane sapphire sub- 
strate, as shown in Fig. [T] We fabricated the junctions by 
e-beam lithography and double angle evaporation of alu- 
minum using the bridge- free technique of |18l 119] . Prior 
to aluminum deposition, the substrate is cleaned of re- 
sist residues using an oxygen plasma |20j . We minimize 
the width of the connecting wires between junctions in 
order to reduce parasitic capacitances to ground, which 
ultimately lower the self-resonant frequency of the su- 
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FIG. 1. (a) Schematic representation of an array of Joseph- 
son junctions. Capacitances of islands to ground and across 
tunnel junctions are denoted with Co and Cj, respectively. 
The junction inductance is given by Lj. (b) SEM image of 
the Josephson junction array fabricated using the bridge-free 
technique 18, 19 . (c) Optical image of an LC resonator. The 
large pads implement the resonator capacitance as well as cou- 
pling capacitances to the CPW feedline. An array of Joseph- 
son junctions between the pads implements the superinduc- 
tance. The ground plane is patterned with flux-trapping 
holes, (d) Low- frequency model for the device shown in (c). 
The phase slip element (split diamond) in series with the su- 
perinductance represents the collective contribution of phase 
slips through all junctions in the array. The characteristic 
energy of the phase-slip element is denoted by Es- 



We characterize our superinductances at low temper- 
atures and microwave frequencies by incorporating them 
in lumped element LC resonators capacitively coupled to 
a co-planar waveguide in the hanger geometry. The res- 
onator response is measured in transmission. An optical 
image of a typical device and its low-frequency circuit 
model are shown in Figs. [ljc,d). 

The samples were mounted on the mixing chamber 
stage (15 mK) of a dilution refrigerator inside a cop- 
per box, enclosed in an aluminum-Cryoperm-aluminum 
shield with a Cryoperm cap. We used two 4-12 GHz 
Pamtcch isolators and a 12 GHz K&L multi-section low- 
pass filter before the HEMT amplifier, and installed cop- 
per powder filters on the input and output lines. 

The internal quality factor of a resonator is extracted 
by fitting the transmission data about the resonance with 
the response function [32] 
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where Qtot and Q ex t are the total and external quality 
factors, fa is the resonant frequency, and Sf character- 
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FIG. 2. (a) Typical microwave transmission data for an 80- 
j unction resonator measured with order one photon circulat- 
ing power; solid line is theoretical prediction corresponding to 
Qint > 37,000. Inset: Temperature dependence of the inter- 
nal quality factor of the resonator. Dashed line is a guide for 
the eye. (b) Same measurement as in (a) for a 160-junction 
array loop; solid line corresponds to Qint > 56, 000. Param- 
eters for the low-frequency model in Fig. [ijd) are C c = 1.6 
fF, Cr = 7.2 fF and L R = 150 nH for (a) and C c = 1.8 fF, 
Cr = 11 fF and L R = 76 nH for (b). 



izes asymmetry in the transmission response profile. The 
internal quality factor is given by Q int = q^3q^ t ■ We 
show a typical measured response for a resonator with 
an 80-junction superinductance in Fig. [2]ja). This yields 
an internal quality factor of 37,000 for the resonator at 
15 mK (stage temperature), corresponding to a loss in 
the superinductance of better than 27 ppm. We note 
that an unknown portion of the internal loss comes from 
the capacitors. The inset of Fig. [2ja) shows the depen- 
dence of the internal quality factor on the stage temper- 
ature. A similar resonator with two 80-junction arrays in 
parallel was measured to have a quality factor of 56,000 
(Fig. [2jb)), corresponding to a superinductance loss of 
better than 18 ppm. The external quality factors, Q ex t, 
of both resonators were 5,000. 

We next evaluate the self-resonant modes of an N- 
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junction array. The Lagrangian of the array is 
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where 4„ are the node fluxes associated with each su- 
perconducting island. We express each node flux as a 
superposition of discrete Fourier mode amplitudes, 
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which leads to a diagonal Hamiltonian of the following 
form: 

N/2 

n= E c^Q fc g_*+L-* (4) 

k=-N/2 

Here Qk are the canonical conjugate "charge" variables 
to the fluxes 4^, and 
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are respective capacitance and inductance matrices in the 
Fourier basis. This immediately leads to a dispersion 
relation of the form, 
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where = 1/ \/LjCj is the plasma frequency of a single 
junction. The coupling capacitance pads at each end of 
the array (see Fig. [ljc)) load down the eigenfrequencies 
(see supplementary information for calculation details). 

The data presented in Fig. [2] corresponds to measure- 
ments of the k = 1 mode. The frequencies of higher 
modes, k > 2, of the array lie outside the band of our 
measurement setup. In order to observe these modes we 
exploit their cross-Kerr interaction, which is induced by 
the junction nonlinearity. This interaction leads to a fre- 
quency shift of the lowest mode (k = 1) when higher 
modes are excited. In Fig. [3^ a) we show the results of a 
two-tone measurement of an 80-junction resonator (same 
device as presented in Fig. |2ja)), where we continuously 
monitor the k = 1 mode while sweeping the frequency of 
a probe tone. When the probe tone is resonant with a 
higher mode of the array, we observe a drop in the k = 1 
frequency. 

Fig. [3jb) shows the comparison between the measured 
frequencies of the array modes and the theoretically pre- 
dicted values. The black crosses show the renormalized 
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FIG. 3. (a) Shift of the lowest frequency mode of the 80- 
junction resonator upon application of a probe tone. Due 
to the weak non-linearity of the array (theory predicts of or- 
der 10 MHz/photon |23|). shifts occur when the probe tone 
matches a plasmonic resonance of the array (marked with red 
arrows), (b) Measured plasma mode frequencies of the 80- 
junction array: The blue square represents the fundamental 
resonant frequency of the resonator (4.355 GHz), and red cir- 
cles indicate the plasma modes pointed out by red arrows 
in (a). The black crosses denote the calculated plasma fre- 
quencies. The gray curve represents the dispersion relation 
without including corrections due to the coupling capacitors 
(open boundary conditions), and markers show corresponding 
plasma frequencies. The inset in (b) shows the voltage profile 
along the array for the first three plasma modes. 



mode frequencies calculated after incorporation of the 
corrections due to coupling capacitance pads (see supple- 
mentary information) , which are in good agreement with 
the measured frequencies represented by colored markers. 
The parameters extracted from the fit were Cq = 0.04 
fF, Cj = 40 fF and Lj = 1.9 nH, with a confidence 
range of 20%. The simulated value of the capacitance 
to ground Cn = 0.09 fF agrees within a factor of 2 with 
the inferred value from the fit. Room temperature resis- 
tance measurements of the junction arrays yield a value 
of Lj = 2.1 nH, which agrees with the fit within 10%. 
Using the fit parameters wc calculate the dispersion rela- 
tion for the bare superinductance, shown as a gray curve 
in Fig. |3jb) . We note that the lowest resonant frequency 
of the bare superinductance is 14.2 GHz, which meets the 
design criterion of having self resonances well above the 
frequency range of interest (1-10 GHz). 

In order to characterize the phase slip rate of the su- 
perinductances, we monitor the frequency of the k — 1 
mode of a resonator formed by two supcrinductances in 
parallel, while sweeping an external magnetic field. As 
flux bias is increased, the persistent current induced in 
the loop increases, and results in a drop in frequency of 
the mode (Fig. [4](a)) . This behavior can be described 
to lowest order in junction nonlinearity by the quasi- 
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FIG. 4. (a) Lowest mode frequency versus applied flux bias 
for resonator with 160-junction array loop. Flux bias is swept 
up (red) and down (blue) over the course of several hours. 
Gray curves represent quasi-classical predictions for resonator 
frequencies at different integer values of flux quanta in the 
array loop. The only adjustable parameter is the number of 
junctions, found to be 150 ± 10. (b) Same measurement as in 
(a) over a smaller flux range, but before each measurement 
a high power pulse at one of the plasma modes is applied in 
order to reset the resonator to the lowest flux state. 
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where <fr e xt is the applied flux bias, and m is the inte- 
ger number of flux quanta inside the loop. A phase slip 
event is associated with an integer change in m, which 
we observe as a jump in resonator frequency. As we in- 
crease the flux bias, the phase slip rate is enhanced and 
frequency jumps become more probable [23]. We swept 
the flux bias applied to the loop over several flux quanta 
before a phase slip event occurred. The typical dura- 
tion between phase slips recorded in this experiment was 
over an hour, Fig. |4]ja). This is a remarkably low phase 
slip rate of well under 1 mHz for a loop of 160 junc- 
tions, significantly lower than previously reported values 
on shorter arrays [H [25] . 

Due to the extremely low phase slip rate, in order to 
measure the resonator in the lowest flux state, we employ 
an active resetting scheme. We apply a high power pulse 
at one of the k > 2 modes before measuring the loca- 
tion of the lowest resonant frequency. This high power 
pulse activates phase slips, allowing the loop to settle 
into the lowest flux state. Using this protocol, we tracked 
the resonant frequency in the lowest flux state, as shown 
in Fig. [4jb). Discrete inverted parabolas are observed, 
which allow us to unambiguously calibrate the number 
of flux quanta in the loop. The fitted value for the to- 
tal number of junctions differs by 6% from the actual 
number, which can be explained by the classical nature 
of the theory which does not take quantum fluctuations 
into account. 

In conclusion, we have demonstrated the successful im- 
plementation of superinductances using arrays of Joseph- 
son junctions, and performed the first detailed charac- 
terization of a superinductance. We measured superin- 



ductances in the range of 100-300 nH, self resonant fre- 
quencies above 10 GHz, internal losses less than 20 ppm, 
and phase slip rates below 1 mHz. The long lifetimes of 
persistent current states in the array loop also demon- 
strate the low DC resistance of the arrays. With these 
parameters, the Josephson junction array superinduc- 
tance enriches significantly the quantum electronics tool- 
box. Applications would include further suppression of 
offset charges in superconducting qubits [5], high Q and 
tunable lumped-element resonators, on-chip bias tees, ki- 
netic inductance particle detectors [26] , and precise mea- 
surements of Bloch oscillations [37J . 
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implementation of a superinductance using a different ar- 
ray topology [2"5] . 
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SUPPLEMENTAL MATERIAL FOR 
IMPLEMENTATION OF SUPERINDUCTANCES FOR QUANTUM CIRCUITS 

SIMULATION OF PARASITIC CAPACITANCE IN JUNCTION ARRAYS 



The parasitic capacitance to ground Co of the islands in a Josephson junction (J J) array lowers the self-resonant 
frequency of the supcrinductance. We optimize the array geometry in Ansoft Maxwell electromagnetic field simulation 
software in order to minimize the parasitic capacitance. 

Our model for the array includes the connection pads and substrate (Fig. |Sl[a)). We model the junction oxide 
layer by a 1 nm thick dielectric material with relative permittivity of 6. Despite only simulating 5 junctions in the 
array, we can robustly obtain the parasitic capacitance for the full structure. Increasing the number of junctions in 
the model leads to no significant change in the computed value of the parasitic capacitance. 

Numerical simulations revealed that junctions with large, 1:30, aspect ratios minimize the parasitic capacitance 
(Fig. |Sl[b)). Fabricating junctions with such large aspect ratios using the Dolan bridge technique is limited by the 
collapse of the bridge. The bridge-free technique (BFT) jTSJ, Q15] allows the fabrication of arbitrary size and aspect 
ratio junctions. The BFT also minimizes the width of the connecting wires between junctions. Our simulations show 
that this leads to a reduction of 60% in parasitic capacitance as compared with the standard Dolan bridge technique. 

The parasitic capacitance Co of each island limits the maximum number of junctions in the array to 



N„ 




(SI) 



a value such that the first self-resonant mode of the array lies ~30% below the plasma frequency of a single junction 
[25] , In Fig. SI b) we present a plot of the simulated Co and N max as a function of JJ size for large aspect ratio 
junctions. For our device we obtain N max = 70. 
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FIG. SI. (a) Simplified BFT junction array model. The orange color represents the thin layer of dielectric (1 nm thick 
AlOa; oxide), and purple represents the perfect conductor (superconducting aluminum), (b) Simulated parasitic capacitance 
to ground and inferred N max (Eq. (Sll) for the design in (a). The junction area is varied by sweeping its length. The dashed 



vertical gray line indicates the parameters for the sample presented in the main text. The lines joining points are a guide for 
the eye. 



CALCULATING THE SCATTERING MATRIX FOR THE HANGER GEOMETRY 



In this section we describe the calculation of scattering matrix of an array resonator coupled to a CPW in a hanger 



geometry. The circuit model is shown in Fig. S2^a), which corresponds to the device shown in Fig. SI c). Impedance 
Z\ represents the coupling between the CPW feedline and the nearest pad. Impedance Z^ represents the cross 
capacitance between the feedline and the opposite pad, while Z 5 and Z 3 represent coupling of the pads to ground. 
Impedance Z^ represents the JJ array. 



The circuit model in Fig. |S2[a) can be simplified by applying the standard Y-A transform, as shown in Fig. S2 c 
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FIG. S2. (a) Model of resonator coupled to a CPW in the hanger geometry. The schematic in (b) elucidates the Wheatstone 
bridge topology of (a). The resultant circuit after a Y-A transform of the circuit in (b) is shown in (c). 



The resulting impedances Z[, Z' 2 , Z' 5 are given by: 

yj Zi Z\ +Z\ Zi (S2) 

yj Z-2 '/.-,-• '/.x> Z\ +Z\ Z2 (S3) 

z > 5 = z^+zj.+z.z^ (S4) 

The total impedance of the structure is given by 

Z 1 Z 2 Z 3 + Z X {Z 2 + Z 3 )Z A + {{Z x + Z 2 )Z 3 + (Zj +Z 2 + Z 3 )Z 4 )Z 5 



^total 



Z 4 (Z 3 + Z 5 ) + Z X {Z 2 +Z 3 + Z 5 ) + Z 2 (Z 3 + Z 4 + Z 5 ) 



(S5) 



We note that the frequency response of Z tota \ near resonance can be approximated using a simplified model as that 
shown in Fig. 1(d) of the main text. Given the total impedance Ztotai which shunts the CPW feedline, the scattering 
matrix can be directly computed using standard techniques |30j . 



DETAILS ON FABRICATION PROCEDURES 



The JJ arrays are fabricated by e-beam lithography at 100 kV using a Vistec 5000+ electron beam pattern generator. 
The use of a high energy electron beam minimizes forward scattering in the PMMA/MMA resist bilayer and allows 
implementation of a bridge- free double angle evaporation technique [18] . We develop using a mixture 1:3 of IPA:water 
at 6 °C. Prior to aluminum deposition, the substrate is cleaned of resist residues using a high pressure, low power 
oxygen plasma 20J. The aluminum films are deposited in a Plassys UMS300 UHV multichamber e-beam evaporation 
system with a base pressure of 10 -9 Torr. The static oxidation in between the aluminum layer depositions is performed 
in a separate chamber, in a mixture 1:3 of O^Ar at 100 Torr for 10 minutes. We obtain critical current densities of 
the order of tens of amperes per square centimeter. The junctions age by <10% during the first week after fabrication. 



MODIFICATION OF ARRAY MODE FREQUENCIES BY COUPLING PAD CAPACITANCES 



The end capacitance pads lower the eigenmodes of the JJ array. To evaluate the modified mode frequencies of the 
./V-junction array k € [—N/2, N/2], we model each k mode as a transmission line of characteristic impedance, 
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FIG. S3. Transmission line model for k modes of an JV-junction array shunted by pad capacitances Cs- 



terminated by a shunt capacitance C s at both ends. Here R q = h/(2e 2 ) = 1.02 kil represents the reduced supercon- 
ducting resistance quantum (Fig. S3 ) . 
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FIG. S4. Determination of fk=i and /fc = 2 modified by pad capacitances, using Eqs.(S8l and |S9| respectively. The fitting 
parameters were C /Cj = 0.001, C s = 0.7 fF and Lj = 1.9 nH. 



In order to calculate a mode frequency of the loaded transmission line, we exploit the symmetry of the problem by 
considering a half-section of the line. The input impedance of each section can be written as |30| 



Z, 



Z L + jZ k tan(/3 fc JV/2) 



(S7) 



Z k +jZ L ttm(J3 k N/2y 

where Zl = l/juJ k C s and /3k denotes the wave number associated with the k th mode. At resonance frequency oj£, 
the impedances of left and right sections of the transmission line satisfy the identity Zj° lt — [Z^ ght ]* . Further, due to 
symmetry, the input impedances of the left and right sections of the transmission line are equal, Z\° lt = Zj" ght . These 
two conditions imply Z- m = 0, which leads to a relationship, 
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Here v® — Nuj®/(kir) denotes the phase velocity for each mode evaluated using the unloaded array calculation. 
Equation (S8) is a transcendental equation in u> k , which can be solved numerically for each of the k modes to obtain 



the corrected mode frequency. The above calculation yields the corrected mode frequencies cj£ only for the odd modes 
of the array. Similar treatment using the input admittance Yi n yields the corrected frequencies for the even modes, 
as solutions of the following transcendental equation 



u%C s Z k = tan( ^ — 



(S9) 



Fig. S4 shows example calculations for the first two mode frequencies using the method described above. 
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